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^ , Abstract 

■ Integer moments of the spectral determinant | det{zl — W)]"^ of complex random 

, matrices W are obtained in terms of the characteristic polynomial of the Hermitian 

I matrix WW* for the class of matrices W = AU where yl is a given matrix and 

J> ' U is random unitary. This work is motivated by studies of complex eigenvalues of 

. random matrices and potential applications of the obtained results in this context 

Q I are discussed. 

(N ■ 

o : 

Q ; 1 Introduction 

^ ; 

P_|, Characteristic polynomials of random matrices have recently attracted considerable in- 
fH I terest in the mathematical physics literature. Initially, the interest was stimulated by 
I applications in number theory [331 EH], quantum chaos O |2I1 HZI and quantum chro- 

^ ■ modynamics (QCD) |121 1201 122], but with the emerging connections to integrable 
^ ■ systems [SHI EE], combinatorics [inillSl, representation theory 011011121113 and analysis 

• th , it has become apparent that characteristic polynomials of random matrices are also 

rS I of independent interest. 

d I In this paper we are concerned with integer moments of the squared modulus of char- 

acteristic polynomial of complex random matrices in a rather general class of matrices 
W = AU, where A a given positive semidefinite Hermitian matrix, A > 0, and f/ is a 
random unitary matrix distributed uniformly over the unitary group. 

In the particular case when A is identity matrix, the matrix W is random unitary, and 
its eigenvalues lie on the unit circle. Various moments of the characteristic polynomial 
for this class of matrices were obtained recently, see jSSl IHHl UHl HH UH] • In the general 
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case, the eigenvalues of W = AU will be distributed in a region in the complex plane. 
Eigenvalue statistics of such complex eigenvalues, and in particular the mean eigenvalue 
density, are of interest for physics of open chaotic systems, see, e.g. 1211122], and in QCD, 
see, e.g. ^] and references therein, and are difficult to study analytically. In this context 
moments of the squared modulus of the characteristic polynomial frequently provide a 
very useful tool. Indeed, in a variety of random matrix ensembles the mean eigenvalue 
density, 

p{x, y) = {ti S{zl -W))w, z = x + iy, (1.1) 

can be expressed in terms of the mean-square-modulus of the characteristic polynomial in 
a closely related random matrix ensemble. In ()1.1|) the angle brackets stand for averaging 
over the matrix distribution, and / is identity matrix. 

An obvious example is served by the Ginibre ensemble of complex matrices jSIj. In 
this ensemble the matrix distribution has density Const, x exp{—tTWW*) where W* is 
complex conjugate transpose of W. The mean density pn{x,y) of eigenvalues of Ginibre 
matrices of size n x n is given by 

1 ^ |^|2fc 

J- _U|2 \Z 



vr k\ 

k=0 

One can arrive at ()1.2|) in various ways. Ginibre computed the joint probability density 
function of all eigenvalues and then applied the method of orthogonal polynomials. An- 
other way is to use the method of dimensional reduction, see e.g. jUlIlHITHj, which yields 
the following relation 

e-l^l' 2 

Pn{x,y) = — — (|det(z/„_i - Wn-i)\ )w„^i- (1-3) 

7r[n — 1)1 

Here the angle brackets stand for averaging over the Ginibre ensemble of complex matrices 
of size (n — 1) x [n — 1). The mean-square on the r.h.s. in ()1.3|) can be easily computed 
yielding again ()1.2|) . 

A less obvious example, which in fact provided the initial impetus for the present 
study, is the so-called ensemble of 'random contractions' j2SI- In its simplest variant of 
rank-one deviations from the unitarity, these are random n x n matrices satisfying the 
constraint 

WnW:=(^^~^ 0<7<1. (1.4) 

In the 'polar' coordinates, Wn = GnUn where f/„ is a CUE„ matrix, i.e. a matrix drawn at 
random from the unitary group U (n), and Gn = diag (a/1 — 7, 1, . . . , 1). The mean density 
of eigenvalues^ of Wn can be expressed as the mean square modulus of the characteristic 
polynomial of (n — 1) x {n — 1) matrices G^-iUn-i, 

p„(x, y) = (^] (I det(z/„_i - Gn-iUn-i)\')u„_„ 1 - 7 < kl' < 1, (1-5) 

7r7|z|^ V 7 J 



^Note that constraint (jl.4|l implies that the eigenvalues of W lie in the annulus 1 — 7<|zp<l. 
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where now the angle brackets stand for averaging over the unitary group U{n — 1) with 
respect to the normahzed Haar measure, and 



Gn-i = diag(A/l - 7, 1, . . . , 1), 



7 



\z\ 






\z\ 
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Another example is provided by finite-rank deviations from Hermiticity jzo] . We only 
consider the simplest but still non-trivial case of rank-one deviations. Let 

Wn = Hn + lTn (1.6) 

where is a GUE„ matrix, i.e. random Hermitian matrix of size n x n with probability 
distribution 



dPp,n{H) = Const. X e-^'^''" l[dHj, J] dH,,dH,, ^ ^ ^ 

j=l l<j<k<n 

and r„ = diag(7, 0, . . . , 0), 7 > 0. It is apparent that all eigenvalues of W lie in the strip 
< y < 7. For the mean eigenvalue density p„(x, y) of Wn we have 

pn{x,y) =rp^n{x,y) {\det{zln-i- {Hn-i + iTn-i)\^)H,,_-,, 0<?/<7, (1.7) 

where 

r,,„(x, y) = 4^^._i(,_2), ' ^.-i = diag(7 - 0, . . . , 0), 

and the angle brackets stand for averaging with respect to the distribution dPp^n-i{H). 
We derive (jl.Sp and ()1.7p in Section IHl 

The above formulas relating the mean eigenvalue density and the mean-square-modulus 
of characteristic polynomial are specific to the considered matrix distributions. In general 
situation, the mean density of eigenvalues can be determined from fractional moments 
of the squared modulus | det{zl — W)\'^ = det{zl — W){zl — W)* (e.g., by the way of 
the logarithmic potential of the eigenvalue distribution), or from averages of ratios of 
det[{zl - W){zl - Wy + e^/], see e.g. plj. Getting explicit formulas for that kind 
of objects outside the classes of Hermitian and unitary matrices is, however, a consid- 
erable challenge. Although it is known that the fractional moments of | det(2;J — W)\'^ 
can be written in terms of a hypergeometric function of matrix argument WW* the 
corresponding series is hard to deal with in the limit of infinite matrix dimension. 

Our main result. Theorem^ expresses 



([det(z/ - AU){zI - AUyf'^)u, m = 1, 2, 



where the integration is over unitary matrices U with respect to the Haar measure, as an 
m-fold integral of powers of the characteristic polynomial of the (Hermitian!) matrix AA*. 
This integral can be written an m x m determinant with entries given by a certain integral 
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transform of the characteristic polynomial of AA*, see ()2.11|) - ()2.12|) . In particular, this 
result implies that for the ensembles of random complex matrices W with unitary invariant 
matrix distribution (for example, for the Feinberg-Zee ensemble ^H] whose probability 
density of matrix entries of W depends only on WW* ) our formulas effectively reduce 
the original non-Hermitian problem to a Hermitian one, albeit on the level of characteristic 
polynomials. This, as explained in more detail at the end of the next section, has a clear 
computational advantage, as one can then use various formulas for averages of products 
and ratios of characteristic polynomials of Hermitian matrices which have been obtained 
in recent years, see jTH EEl ISj- In contrast, with the exception of essentially Gaussian 
weights [HE], no such formulas are known for complex matrices. 
We also express 

I \ 

det[{zl - AU){zI - AUy + eH] / ^ 

as a two- fold integral of the inverse spectral determinant of AA*, see Theorem |21 Again, 
the non-Hermitian problem is reduced to a Hermitian one. This regularized inverse spec- 
tral determinant can be useful as an indicator of the domain of the distribution of complex 
eigenvalues. For, any point z where the above average blows up in the limit £ — > must 
belong to the domain of eigenvalue distribution. 

Acknowledgments. We would like to thank A. Gamburd for useful discussions and, 
in particular, for bringing reference |13] to our attention and for pointing us towards the 
link between our Lemma El and the Selberg Integral. We are also grateful to Ph. Biane 
for bringing reference jlHj to our attention. 



2 Statement of main results and discussion 

Let n and m be positive integers. Define 

^ m 

dflnih, . . . , t™) = - A2(ti, . . . , t^) + tj.) — 2™ dh... dt^, tj > 0, (2.1) 

and, for n > 2m, 

^ m 

dUn{tl, ...,tm) = Y A\t,, ...,tm) ' t,)""'"^ d^ . . . dt^, < tj < 1, (2.2) 

where 

A(ti,...,U = det(trOr.=i= n (^^-^^■) (2.3) 

l<i<j<m 

is the Vardermonde determinant, and 

n + m+? ! m — ? — 1 ! 
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are the normalization constants. The Selberg Integral, see e.g. jHE], asserts that dfin and 



dun are unit mass measures, 



oo oo 



dUnih, . . . ,tm) = j . . .j dUnih, . . . , t^) = 1. 


The measures dfi^ and dun define probability distributions which have the following 
random matrix interpretation. Consider two families of matrix distributions on the space 
of m X m complex matrices Z = {xjk + WjkjJ'k^i- 

1 1 

dfiniZ) = ^ det"+'"(/^ + ZZ*) n ^^^^ ^^J-fc' ^ ^ 0' (2.5) 

j,k=l 



and 

dUriiZ) = 

k. 



dPn{.Z) = i det"-2"(/^ - ZZ*) f[ dxjk dyjk, n > 2m. (2.6) 

j,k=l 

The measures dOniZ) are defined on the matrix ball ZZ* < 1^ and the constants c„ and 
kn are determined by the normalization condition 



/ djln{Z) = / duniZ) = 1. 

JZZ*>0 JQ<ZZ*<Irr, 



A standard calculation, see e.g. jSOj, shows that if s{ZZ*) is a symmetric function of the 
eigenvalues ti, . . . ,tm of ZZ*, i.e. s{ZZ*) = s(ti, . . . , tm), then 

+ 00 +00 

S{ZZ*) dflniZ) = I ■ I S{ti,...,tm) dfin{tl,...tm), (2.7) 

ZZ*>0 J J 



1 1 

[ s{ZZ*)dVn{Z) = [■■■ [ S{h,...,t^) dUniti,...tm). (2.8) 

J0<ZZ*<Im J J 

~ ~ 

Theorem^below, which we state in a slightly more generality than required for spectral 
determinants, tells how to integrate moments of determinants over the unitary group 
equipped with the Haar measure dU fixed by the normalization /^^^^ dU = 1. 



Theorem 1 Let A,B,C,D be complex matrices of size n x n. 
(i) For any positive integer m 



dei'^[{AU+C){BU+D)*]dU = ... JJ det(CD*+tjA5*) c//i„(ti, . . . , t^) (2.9) 

U{n) ■5''=^ 
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(a) If AA* < CC* and BB* < DD* then for any positive integer m such that 2m < n 

(2.10) 



1 1 

dU f f dUn{tl, . . . ,tr, 



Uin) . 



Remark. Identities ()2.9p - ()2.10j) may be written in yet another form by making use 
of the well-known identity 



det {pj{ti)). .^^ det {qj{ti)). .^-^ dti . . . dt^ = m\ det ( / pi{t)qj(t) dt 



We have 



+00 \ '"-l 



U{n) \0 / ij=o 

and 



\ \ m—l 



dU ml I r (1 - tY-^"" dt , 

det / \„,J^, . (2.12) 



det"" [{AU + C){BU + D)*] fc„ \J det{CD* - tAB*) 

U(n) \0 / ij=o 

Obviously, by letting C = D = 2;/ in ()2.9|) and ()2.im) one obtains formulas for moments 
of the spectral determinants | det(2;/ — AU)\'^. In particular, 



00 

2 



. M7 / X , deti\z\^I + tAA*) , 
det(2/ - AU)\' dU ={n + l) / J^^^ dt (2.13) 



U{n) 

and, provided n > 2, 



dU 



l + t) 



( 1 

n-2 



C/(n) 



det(z/- Af/)|2 





/■ (^ - 1)(1 -t)""^ , , ,0 , , , 



(2.14) 



where Amm(^^*) and Amaa;(AA*) are respectively the smallest and largest eigenvalues of 
AA*. 

If \min{AA*) < |zp < X-maxiAA*), then the integral on the left-hand side in p.l4j) 
should be handled with care. One way to do this is to regularize the integrand. 
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For positive e, define 

dU 



n) det [eH+{l-\AU) {l - \AUy] 



The integral on the right-hand side is, in fact, a function of A A* and our next theorem 
evaluates this function in terms of the eigenvalues of AA* . 

Theorem 2 Let e > 0, and assume that n > 2. Then for any n x n matrix A and any 
non-zero complex z 

n-1 f\_ f^'^dx 1 



i?,,,(A,A*) = — - / {l-tr-'dt 

Jo J~oo ^ det 



^AA* + {e^-t)I- leVt [x + I) I 



If the eigenvalues a| of AA* are all distinct, then for any non-zero z in the annulus 
\^in{AA*) < < X„,ax{AA*) we have 

^ ' ' j=l k^j ^ i 

where 9 is the Heaviside step function, 6{r) = 1/2 at r = 0. 

We prove Theorems Q and El in Sections 0] and respectively, by making use of two 
techniques, which to a certain extent are equivalent. One is based on the expansion 
of moments of spectral determinants in characters of the unitary group and subsequent 
use of the orthogonality of characters. On this way. Theorem Q is equivalent to two 
combinatorial identities (jH.lSj) and one of which is a particular case of the Selberg 

integral in the form of Kaneko jM| and Kadell ^S] • We prove ()H.18|1 and in Section 

EJ These combinatorial identities can be stated in the form of matrix integrals ()3.28|) 
and ()3.29|) and are of independent interest. They lead to evaluation of some non-trivial 
matrix integrals, as discussed at the end of Section El The other technique is based on 
the so-called color- flavor transformation, due to Zirnbauer |1H|- This transformation has 
many uses, and in the random-matrix context it provides a very convenient tool to handle 
moments of spectral determinants. 

As an application of Theorem ^ let us consider random matrices ()1.4|) . In the limit 
n ^ oo the eigenvalues of Wn get closer and closer to the unit circle. Let Mn{a, b) be the 
number of eigenvalues of Wn in the annulus 

.2a , ,p 2b\ 

D„h= < z : — < 1- z(^ < — y 0<a<b. 



n n 
By (USD, 

{Mn{a,b))ui^n)= j ^^(^) / \^^^{zIn-l~Gn-lUn-l)?dUn^l 

Da,b \ U{n-1) J 



dxdy. 
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Making use of (jTT!?|l . 



U{n-1) 

and 
where 



1 

-{J^nia,b))u{n) = TT / fn{,q)dq, 

n Ji_2b 

n 

' »-7« Jo (! + ()"+' 

and 7 = (|zp + 7 — Letting n — oo, we obtain, after simple manipulations, that 

1,.,. sinha /a(7-2)\ sinhfc fb{-f-2)\ 

hm -(A/„(a,6))[/(„) = exp ^xp , 

n^oo n a \ 1 J b \ 7/ 

recovering one of the formulas of j2SI, who, using a different method requiring knowledge 
of the joint probability distribution of eigenvalues, found the mean density of eigenvalues 
and higher order correlation functions for the general case of finite-rank deviation from 
the CUE. Note that when 7 = 1 the nonzero eigenvalues of GnUn coincide with the 
(n — 1) X (n — 1) matrix obtained from f/„ by removing its first row and column, see [IHj 
for more information about eigenvalue statistics of truncated unitary matrices. 

Now, we would like to elaborate on the point made at the end of Introduction. Consider 
random complex matrices W of the size n x n with unitary invariant matrix distribution. 
Then, by making use of the unitary invariance and Theorem ^ 

{\det{zI-W)\')^ = ( [ \det{zI-WU)\'du) = H .^"i'^'l^ dt, (2.16) 

\Ju{n) I ^ Jo U + 

where 

Vn[x) = {det{xl + WW*))w. 

A similar formula holds for higher order moments of | det(2;/ — Thus, Theorem^ 

reduces the original non-Hermitian problem to a Hermitian one. 

The integral on the right-hand side in ()2.1(ij) can be evaluated, in the limit of infinite 
matrix dimension, in terms of the limiting eigenvalue distribution of the non-negative 
matrices WW*. To this end, consider, for example, the complex n x n matrices W with 
the matrix distribution characterized by the Feinberg-Zee density 

Const. X e-«t'^^(^^*), (2.17) 
where V{r) is a polynomial in r, V{r) = amf"^ + . . ., > 0. Then 

p„(x) = e"^'"("+^)''"'(^)(l + o(l)). 
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where dw{X) is the hmiting normahsed eigenvalue counting measure of WW*, and it can 
be shown that 

hm -ln(|det(zJ-W^)p) =$(x,i/), 



n^oo n 

where 



In if \z\ > nil = J Xdw{X), 

\nXdw{X) if ^ > m„i = /" (2.I8) 



z\'^ + I \n-—-^ — —dw{X) if l/m_i < |z| < mi 



12 , / i„ ^ + ^0 

'0 kP + ^o 



where to is the unique non-negative solution of 

dw{X) 1 



X + t \z\'^ + t 

The function $(a;,y) is subharmonic and, hence, defines a measure du = in the 

complex plane. Here A is the Laplacian in variables x and y. For the Ginibre ensemble 
of random matrices this measure can be found explicitly. In this case V{r) = r and 
WW* is a Wishart ensemble of random matrices. Its limiting eigenvalue distribution 
dw{X) is given by dw{X) = {A — X)/ X, < A < 4, with mi = 1 and m_i = 00. A 
straightforward but tedious calculation shows that ^{x,y) = — 1 inside the unit disk 
< 1. Therefore du is the uniform distribution on the unit disk, which is the same 
as the limiting eigenvalue distribution in the Ginibre ensemble of random matrices, and, 
hence, for this ensemble 



Jim i (In I det(./ - W)\X = i In ( | det(./ - W)\X , (2.19) 

so that the operations of taking logarithm and taking average commute in the limit 
n ^ 00. A similar relation is known to hold for Wigner ensembles of Hermitian matrices 

m 

It would be interesting to investigate conditions on random matrix distributions which 
guarantee ()2.19j) . As the left hand-side in ()2.19|1 is the logarithmic potential of the limiting 
eigenvalue distribution of W, this together with our Theorem [T] would give a useful tool 
for calculating eigenvalue distributions in the complex plane. There are indications that 
the range of matrix distributions for which ()2.19|) holds is quite wide and contains the 
invariant ensembles ()2.17|1 . Indeed, ^{x,y) of equation ()2.18|1 reproduces the density of 
eigenvalue distribution in ensembles ()2.17|1 which was obtained in fOl 120] with the help 
of the method of Hermitization^. In this context we would like to mention calculation of 
Brown's measure for i?-diagonal elements in finite von Neumann algebras [28], see also jH]. 

^This method has a hidden regularization procedure which has to justified to satisfy the mathematical 
rigor. 
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A matrix model for such elements is provided by random matrices RU where U is random 
unitary and R is positive-definite, and Brown's measure is in a way a regularized version of 
the eigenvalue distribution. Again, y) of equation ()2.18|) reproduces Brown's measure 
found in pSj. 

3 Combinatorial identities 

Schur functions. In order to make our paper self-contained we recall below the required 
facts from the theory of symmetric polynomials. 

A partition is a finite sequence A = (Ai, A2, . . . , A„) of integers, called parts, such that 
Ai > A2 > . . . > An > 0. The weight of a partition, |A|, is the sum of its parts |A| = J^j ^j^ 
and the length, /(A), is the number of its non-zero parts. No distinction is made between 
partitions which differ merely by the number of zero parts, and different partitions of 
weight r represent different ways to write r as a sum of natural numbers. 

Partitions can be viewed as Young diagrams. The Young diagram of A is a rectangular 
array of boxes (or dots), with Xj boxes in the jth row, the rows being lined up on the left. 
By transposing the diagram of A (i.e. interchanging the rows and columns) one obtains 
another partition. This partition is called conjugate of A and denoted by A'. For example 
the conjugate of the partition (r) of length one is the partition (!,...,!) = (1^) of length 
r. Obviously, /(A') = Ai and |A| = |A'|. 

For any partition A of length 1{X) < n 

det 

SxiXi, . . . ,Xn) = : n (3-1) 

det {xr%=, 

is a symmetric polynomial in xi, . . . ,Xn, homogeneous of degree |A|. These polynomials 
are known as the Schur functions. By convention, Sx{xi, . . . , x„) = if Z(A) > n. This 
convention is in agreement with the apparent identities 

Sa(Xi, . . . ,X„_i,0) = SxiXi,...,Xn-l) if /(A) < 72 - 1 (3.2) 

= if/(A)>n-l. (3.3) 

For partitions of length one, A = (r), the Schur functions sx are the complete sym- 
metric functions hr, 

S(r) (3^1 ) • • • ) Xn) hj. (xi , . . . , Xji^ ^ ^ 3^il ^Jjj • • • X^^ , (3.4) 

l<jl<j2<---<jr<" 

and sx' are the elementary symmetric functions e^, 

S{l^){Xi, . . . , Xn) = Gr{Xi, . . . , Xn) = ^ ^ Xi^^Xi^ ■ ■ ■ Xi^ . (3-5) 

l<ii<i2<...<ir<n 
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More generally, see e.g. P- 41, the Jacobi-Trudi identity asserts that for any n > /(A), 



sx = det (/iA,-i+j)^=i , sy = det (eA,_»+j)^^i (3.6) 

where, by convention, Cr = hr = if r < 0. 

We shall also need the Schur functions of matrix argument. If M is an n x n matrix 
then 

Sx{M) = Sx{Xi, ...,Xn) 

where the eigenvalues of M. Thus sx{M) is a symmetric polynomial in the 

eigenvalues of M. In view of fl3.(ij) . it is also a polynomial in the matrix entries of M. The 
Schur functions of matrix argument are the characters of irreducible representations of the 
general linear group and its unitary subgroup and, as a consequence, have an important 
property of orthogonality. If A and /i are two partitions and A and B are two n x n 
matrices then, see e.g. P- 445, 

/ sx{AU)j;XBU)dU = 6x /'^'^^*\ (3.7) 

Ju{n) "A 

and 

/ sxiAUBU*)dU = (3.8) 

Ju{n) "A 

where dx is the dimension of the irreducible representations of U{n) with signature A, 

dx = Sx{In) = Sx{ln)- 

We use the notation (1„) for the n-tuple (1, . . . , 1). 
If A is a partition of length 1, A = (r), then 

n ^ h ^^ \ f n + r-l\ {n + r - 1)\ 
sxiln) = Kiln) = ^ J = , (3.9) 

and 

(T) \ Tl^ 
= u ' M - (3-10) 
r J r\{n-r)\ 

In general, explicit expressions are known for sx{ln) and sy{ln) in terms of the Aj's. If 
/(A) < n then for any m > /(A) 

^a(i„)={ n (A.-^-A,+,)|xn ^^j"^+_\7/j^^.^, . (3.11) 

Kl<i<j<m J 

If Z(A') < n then for any m > /(A) 



L {m + Xj - j)\{n - 



sy{ln) = { n (A.-^-A,+j)l xfj ^ . ^V„/^Lx (3.12) 

'i<i|<^ J ^i(^ + J-l-A,)!(m + A,-j)! 
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Both identities can be derived by evaluating the binomial determinants in ()3.6p . 

We shall also need the Cauchy identities for Schur functions, see, e.g., jSTj, pp. 63, 65. 
Let X be an X n matrix. Then 



fJdet(/„ + t,X) = J2sx{h,...,tJsy{X) (3.13) 

i=l A 

The summation in ()3.13|) is over all partitions such that /(A) < m and /(A') < n and is 
finite. The summation in ()3.14j) is over all partitions such that /(A) < min(m, n) and is 
infinite. The corresponding series converges absolutely if XX* < In- 

Beta-function determinants. When m = 1 identities ()3.13|1 and ()3.14|1 take the 
familiar form of the expansion of the characteristic polynomial and its reciprocal in terms 
of the elementary symmetric functions and complete symmetric functions, respectively. 
In this case Theorem ^ is a straightforward consequence of the orthogonality property of 
the Schur functions ()3.7|) and the Euler integral 

/■+°° fP^^ T('n)T(n) 
/ tP-\l-ty~Ut= dt = -^^-^^ = B(p,q), Rep,g>0, (3.15) 

where T{p) and B{p,q) are the Gamma and Beta functions respectively. Indeed, for 
example. 



„ n 

/ det{I + AU) det{I + BUy dU = V 

Ju(n) 



CriAB*] 
Cr(ln) 



= (n + 1) / det{I + tAB*] 

where we have used ()3.15|) in the form 



dt 



[1 + ty 



+2- 



1 



+ 00 



e.(ln) ^ 'Jo (l + t)"+2 

Similarly, 

1 . 

= (n 



^ + 1)/ T^TZ^dt- (3-16) 



1) / f{l-ty-^dt, (3.17) 
Jo 



and this identity does the trick for the reciprocal characteristic polynomials. 

Our proof of Theorem [T] uses the following generalization of ()3.16p - ()3.17|) to multi- 
variate setting. 



Lemma 3 Let m and n be nonnegative integers. 
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(a) For any partition A such that /(A) < m and 1{X') < n 

^i&4 = l [...[ sx{t,,...t^)A\t^,...,t^)Y\- ^^^f-^. (3.18) 



(b) If 2m < n then for any partition A such that /(A) < m 

11 ™ 



/n II V 

. . . Sxit,, . . . trr,)A\t,, . . . , t„) J](l - (3.19) 

-^^^ 



r/ie normalization constants Cn and kn are given in ^2.4\ ), and A(ti, . . . ,tm) is the Van- 
dermonde determinant \2.^) . 

Remark. Identity ()3.19p can be inferred from a generalization of the Selberg Integral 
due to Kaneko [31] and Kadell of which ()3.19p is a particular case corresponding to a 
special choice of parameters. However, we are not aware about any generalization of the 
Selberg Integral leading to ()3.18j) . Below, we give an elementary proof of ()3.18|1 and ()3.19j) 
which is based on evaluating a determinant consisting of Beta functions, see Proposition 
13] below. Our proof has a limited scope and does not extend to the generality of Kaneko 
and Kadell formulas. 

Proof. Let fj = m + \j — j, j = 1,2,..., m. If /(A) < m < n and 1{X') < n then by 

dHU) - 



where 

A(/i,...,/J= n (/^ - /.) = det (/rX,=i • 

l<i<j<m 

By adding rows in the Vandermonde determinant det 

A(/i, ...,/„) = det {Pm-iUj))Z=i ' 



where 
Hence 
and 



Pk{x) = {x + l){x + 2)...{x + k). 
f.W . . . fJ.Aif,, ...Jm) = det ((/,■ + m - , (3.20) 



'm—l 
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where B is the Beta function. By making use of Proposition 0] below, 

det( B{fj+m-i + l,n+m-fj) )™ = det ( B{fj +m-i + l,n+m- fj )^^.^-^ 

= det 







+ 



n+2m 



It is apparent that 

f+oo \ +00 +00 



,0 / ij=i 

+00 +00 



(it,; 







A (l + t,)-+2-' 



and (ITT^ follows. 

Similarly, if 2m < n and /(A) < m then by (jS.lljl 



and, in view of ()3.20p . 



By Proposition 

det {B{fj +m-i + l,n-2m+i))^^.^^ = det +m-i + l, n-2m+l))^^.^^ (3.22) 

= det( [ rft) ^ (3.23) 

V-'o y i,j=i 

and (pnnil follows. □ 

Proposition 4 For any pi,p2, . . . , Pm o-nd qi,q2, ■ ■ ■ ,qm such that Kepj > m and Re qj > 
— 1 we have 

det ( B{pj - i, qj + i))'^ = det ( B{pj - i, qj + 1) )™ . (3.24) 

Proof. We shall use the identity 

B{p, q-l) + B{p- 1, q) = B{p-l,q- 1) (3.25) 
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and the operation of addition of columns to transform the determinant on the left in 
()3.24|) to the one on the right. 

It is convenient to write determinants by showing their columns. With this convention, 



\B{p - 1, g + 1), B{p - 2, g + 2), . . . , B{p - m + l,g + m-l), B{p -m,q + m)\ (3.26) 

represent the determinant on the l.h.s. in ()3.24|1 . Let us label its columns by numbers 
1, . . . ,m from left to right (so that the leftmost column is column 1). Note a particular 
property of columns in this determinant. As we move from column j to column j + 1 the 
first argument of the Beta function decreases by one 1 the second argument increases by 
1. To be able to refer to this property, we say that columns 1,2, ... ,m are balanced. 

Observing that column 1 has the desired form already, let us perform the following 
operation on columns 2, 3, . . . , m. Starting at column m and working backwards, let us 
add to each column the one that precedes it. In view of ()3.25|1 and the above mentioned 
property of columns, this operation yields 



B{p-1, q+1), B{p-2, q+1), B{p-3, q+2), B{jp-m+l, q+m-2), B{jp-m, q+m-l)\. 



Observing that columns 1 and 2 have the desired form now, and that columns 2, . . . ,m 
remain balanced, we apply our operation again, now on columns 3, . . . , m. This yields 
the determinant 



B(j)-l, q+l), B{p-2, q+1), B{p-3, q+1), B{p-m+l, q+m-3), B{jp-m, q+m-2)\. 



where columns 1, 2, 3 have the desired form and columns 4, . . . , m are balanced. It is clear 
that repeated application of our operation will yield the determinant 



Applications to matrix integrals. In view of integration formulas ()2.7|) and ()2.8p . 
identities ()3.18p and ()3.19|) can be rewritten as: 



where Z are complex mxm matrices, and Im and In are identity matrices of sizes mxm 
and n X n, respectively. The first identity holds for any non-negative integer n and any 
partition A such that /(A') < n. The second one holds for any integer n > 2m and any A. 
Since the above identities become trivial (both sides vanish) for partitions of length > m 
we drop the restriction /(A) < m. 

These two identities lead to several useful matrix integrals. 

Let M be an m X m matrix. Then, for any non-negative integer n. 



\B{p-l,q + l),B{p-2,q + l),B{p-3,q + l),...,B{p~m,q + l)\ 
after the final step. This is exactly the determinant on the right in ()3.25|) . 



□ 




(3.27) 




(3.28) 
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provided /(A') < n, and if n > 2m then 

/ .a(MZZ*) dUZ) = (3.29) 

JZZ*<Irr, Sx{Inj 

for any A. These two integrals follow from ()3.27p and ()3.8|) and the unitary invariance of 
dfin{Z) and dOni^Z). 

If L and M are two m x m matrices then for any non-negative integer n 

f s,iLZ)JjMZ) dfin{Z) = 5x J^^^^*\ (3.30) 

provided /(A') < n and /(/i') < n, and if n > 2m then 

/ sx{LZ)s,{MZ) dO^Z) = 5x /'^^,f*\ (3.31) 

JZZ*<Im Sx{Inj 

These orthogonality relations follow from ()3.7|1 and ()3.28|) - ()3.29|1 . and, in turn, lead to 
Berezin-Hua integrals jHOl IZ| 

det"(/„ + LZ)det"(/„ + MZ)* dfin{Z) = det"(/^ + LM*) 

ZZ*>0 

dl)n{Z) 1 



zz^<i^ det"(/^ - LZ)det"(/^ - MZ)* det"(/^ - LM*] 



n > 2m. 



One only has to recall the Cauchy identites ()3.13j) and ()3.14|1 . 

If P and Q are two n x m matrices and n > 2m then it follows from ()3.7p and ()3.3ip 
that 



sxiPQ*U)s^iPQ*U) dU = / sxiP*PZ)s^{Q*QZ) du^{Z) (3.32) 

U{n) JZZ''<Irr, 

for any A and /i. Identity ()3.32|1 implies that 

I ^tr(PQ*u+u*QPn = I e*'^^^*^^^^*^*^) rfz>„(Z). (3.33) 

Ju{n) JzZ*<Im 

The duality relation ()3.33|) is a particular case of Zirnbauer's color-flavor transformation 
It can be easily obtained from ()3.32|) by making use of the expansion 



e'^^ = Y,cxSx{A). (3.34) 

A 

In fact, ()3.33|) extends to any series g{A) = J2x^^^>^i^)^ 

[ \g{PQ*U)\^ dU= [ g{P*PZ)g{Q*QZ) dv^{Z). 

Ju{n) JzZ*<Im 
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It follows from ()3.34|) and ()3.7|) that the integral over the unitary group on the left- 
hand side in ()3.33p is a function of Q*QP*P. This function can be evaluated explicitly 
in terms of the eigenvalues of Q*QP*P. We would like to demonstrate this in a slightly 
more general setting. 

For square matrices A and B of size n x n define 

F,,{AB*) = [ e*'-(^^+^*^*) dU. (3.35) 

Ju{n) 

If the eigenvalues zl,...,z'^of the matrix AB* are all distinct then |1T] 

F.iAB*) = ^^1^^ X det (.r^/,_,(..)):,., 
where Ik is the modified Bessel function, 

For our purposes, we want to know Fn{AB*) for matrices AB* of low rank, e.g. when 
AB* is rank one. 



Lemma 5 Suppose that AB* has m distinct non-zero eigenvalues zf, . . . ,zf^ and 2m < n. 
Then 

/™ (n-iV \ } } det (g{tiZ^)r. ^ - 



\{n-m- i)\ I j " J /\{zl,. . . ,z;^ 







where g{x) = Iq {2^/x) . In particular, if AB* is rank one and z^ is its non-zero eigenvalue 
then 



F^{AB*) = (n - 1) Jo (1 - t)"-^ dt. 

Proof. It follows from (jO^ and (jSHI) that 

F^{AB*)=J2-^s,{AB* 

The coefficients cx are given by 



(3.36) 



c.=det(-— -- =..(u n 
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see, e.g., jl] and references therein, and 



m— 1 



where as before fj = m + Xj — j. Note that the summation is over all partitions A of 
length < m, or, equivalently, over all /i > /2 > . . . > /m > 0. It follows now from 1)3.211] 
- (IT^ that 

where 

det ftf^)" det (z-^'Y 

yi,''!, • • • trnj — f \2 f \2 

fi>f2>->fm>0 

To complete the proof, recall the following generalization of the Cauchy-Binet formula, 
see e.g. |30j p. 22. If g{x) = Ylif>olf^^ analytic function in the complex x-plane 

then 

det {g{tiXj))"\^^ = 7/, . . . 7/„ det (t('^ ^ det (^x 

/l>/2>...>/m>0 ^ 

By making use of this formula. 



1 

i,j=i \ / i,j=i 



g{tu...,tm) = det (/o (2^/t^))^^^^ , 
and Lemma follows. □ 



4 Proof of Theorem U 

After all the preparatory work of the previous section. Theorem becomes almost evident. 

We first prove ()2.9|) - ()2.10|) for C = D = I. With Lemma El in hand, this becomes a 
routine calculation. Expanding powers of determinants in the Schur functions as in ()3.13|) 
- ()3.14|1 and integrating over the unitary group with the help of ()3.7|1 . one gets 

[ det™[(/ + AU)iI + BUy]dU = ^^7t4 ^^'(^^*) (^-l) 



U{n) 

and 



/det"[(/ - AU){I - BUY] - ^ liZ) ^^-^^ 

U{n) ^ 
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The sum in M.lj) is finite and the sum in ()4.2j) is absolutely converging for any AA* < I 
and BB* < I. Now, by making use of ()3.18|) and ()3.19p . and then ()3.7|) again, one arrives 
at 

det^'iil + AU){I + BUy]dU = . . . Y[det{I + t,AB*) dfin{tu...,t,n) (4.3) 

Uin) 

and 

1 1 

dU f f dunih, . . . ,tm) ^^^^ 



Extending (j43|l and (jOjl to the generality of (jT9|l and (j2.10|l is straightforward. If C 
and D are not degenerate, then 

j def^KAU + C){BU + Dy]dU = det™(CD*) ^ det™[(J + C-^AU){I + D-^BUy]dU 

U(n) U(n) 

and ()2.9j) follows from ()4.3|1 . The assumption that C and D are not degenerate can be 
removed by the continuity argument. □ 



5 Regularization of the inverse determinant 



In this section we employ another approach to the problem of evaluating of negative 
moments of spectral determinants which is to write the determinants as Gaussian integrals 
and then perform the integration over the unitary group with the help of the color-flavor 
transformation. This approach is not new. It was pioneered by Zirnbauer in the context 
of unitary random matrix ensembles. The new element here is that we apply it in the 
general context of complex matrices. 

We shall write the spectral determinant det[(/ — AU){I — AUy] of n x n matrices as 
2n X 2n block determinant 



det[(J - AU){I - AUy] = det[{U* - A){U* - Ay 
and more generally 

det[e'^I + (/ - AU){I - AUy] -- 



i{U* - A) 
i{U* - Ay 



el i{U* - A) 
i{U* - A)* el 
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Proposition 6 Suppose that Re Aj > 0, j = 1, 2. Then for any complex n x n matrix fl 
we have 

"^AI = — - I d^v [ d^w exT) i — \\iv*v + \2W*w + i(w*Q*v + v*flw)]} . 

(5.1) 

The integral on the right-hand side converges absolutely. 

Remark. In this section, we shall use letters in bold face to represent column vectors 
in C". The symbol d'\) will denote the volume element of v in C", 



d'^v = d \ij = Y\ d Re Vj d Im 



Vj. 



Proof. Note that 

\iv*v + \2W*w + i{w*n*v + v*nw) = {v*, w*) ^^^f (^^ ^ 

In view of the singular value decomposition Q = U*ujV, 

All in \ _ fU* \ fXiI iuj \ (U 

in* Aa/y Vo v*)\iuj Aa/y Vo y 

where uo is diagonal matrix of singular values of fi, = diag(c(Ji, . . . ^uOn)., and U and V 
are unitary matrices. Introducing / = Uv and g = Vw, 

. ^ *x All iVL \ ( v\ . *x Ai J iuj\ f f 



Ai ^cuA / fj 
iuj X2J \gj 



Since U and are unitary, d\) = d'^f and d'hv = d'^g. Changing the variables of 
integration in ()5.1|) from v and w to f and g breaks this 2n-fold integral into the product 
of the 2-fold integrals 



2\-l 



- / exp [-Ail/jf - A2|5'jf - i^j{fj9j + fj9j)\ d'^fjd'^gj = (A1A2 + u;-) 

Thus, the integral on the right-hand side in ()5.Hl equals 11^=1 (^1^2 + which is 

obviously same as the determinant on the left-hand side. □ 

Proof of Theorem\^ Obviously, without loss of generality we can put z = 1. Let 
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The integral on the right-hand side converges for any n x n matrix A. 
It follows from Proposition IHl that 

Rs{A,A*) = ^ [d% [ d^w e-[^("*"+^*^)-^(^*^*"+"*^^)l/„(t;*viy*t«), (5.3) 
where, cf. ()3.35p . 

fn{v*VW*w) = I e'^^*U'u+v*U*.v)^jj ^ f tr(W [/+{/* 

Juin) Ju{n) 

By Lemma m 

fn{v*vw*w) = J Jq (2Vtv*vw*u?j danit), (5.4) 

where 

dan{t) = (n- 1)(1 -t)"-2rft 

and Jo is the Bessel function 

i=o 

We have \fn{v*vw*w)\ < 1 for all v and w. This is because |Jo(^)| ^ 1 for all z. 
Therefore we can interchange the order of integrations on replacing /„ in ()5.3|) by the 
integral of ()5.4|1 . This yields 

Re{A,A*) = [ danit) \ f d\} f d\v e-'i.^'^+^''^)+i{^'A''^+^''A^)j^ (2^/tv*vw*w) . 

(5.5) 

In order to perform the integration in variables v and w we shall make use of the integral 
representation 



1 dx 



Jo(2Vpg) = — / — e^(^^+') (5.6) 



which holds any p > and g > and is a particular case of equation 3.871.1 in The 
integral in ()5.(i|l converges because of the oscillations of the exponential function, however 
the convergence is not absolute. 
We have 

Jo (2Vtv*vw*w) =^ r^'^ ^^Vti..*.+^)_ ^5 7) 

V / 2m J_^ X 

Note that by Proposition 



[e — i\/tx)I —A 
-A* 



-1 
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Therefore, on replacing Jq in ()5.5|) by the integral of ()5.7|) and changing the order of 
integrations one arrives at 



(5i 



X 



<\et[AA* + _ t)I - ieVi{x + i)/] 



which is the identity claimed in Theorem |21 It remains to justify reversing the order of 
integrations with respect to x and v, w. Firstly, we will show that the integral on the 
right-hand side in ()5.8j) is well-defined. 

Proposition 7 For any e > and n > 2 the integral in i5.^} converges absolutely (and 
uniformly in A). 



Proof. Let a| be the eigenvalues of AA* so that 



1 1 

X det[AA* + (e2 _ t)I - ie^/t{x + x -LA 



1 

- JJu7(a,-,t) 



3 

where 

w(a,t) 

f 2 

Since \z\ > I Imzl, we have 



+ —t — ie\/t{x + -) 



1 / 

- w[ai,t) 

X 



1 

< 



x^ 



Also, for all < t < £2/2 we have 

1 2 

-2 



"^"^•'^^'^ k2-t + a2| ^i^ 



and for all £ /2 < t < 1 we have 



Therefore the absolute value of the integrand in ()5.8|) is majorated by the function 



1 / 2 ^ 



eVi{l + x^) W 

which is obviously integrable with respect to dcXnit) x dx. □ 

We can now turn to justification of reversing the order of integrations in v, w and x 
in the integral 

J= f'darXt) [ d\, [ d\ju e--(-*-+-*-)+^(-*^*-+-*^-) e*v^(--*-+^). (5.9) 

Jo Jc" Jc" J-oo ^ 
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The corresponding calculation is routine but tedious. First we restrict the x-integration 
to the finite interval 6 < \x\ < 1/6, 6 > 0, reverse the order of integrations, and then 
show that the corresponding tail integrals in are negligible in the limit 6-^0. 
Let 



1 



C" <5<|x|<l/<5 

The absolute value of the integrand is majorated by the integrable function ^^^-<^{'"*'"+'w*'w) 
and therefore we can reverse the order of integrations and then perform the integration 
in v,w. This yields 

Is = dcr„{t) / — . 

Jo Js<\x\<i/5 X det [AA* + (£2 _ t)/ - teVii^x +1)1] 

It follows from this, in view of By Proposition [3 that 

Is= [ daJt) I — ^ F , + oil) 

Jo Aoo ^ det[AA* + {e^-t)I-teVi{x + l)I] 

in the limit 5^0. It only remains to show that the tail integrals 

X/= f'da^it) [ d\ [ d'w e--(-*-+-*-)+^(-*^*-+-*^-) / — e^Vt{x.'v+^) 

Jo Jc" Jc" J\x\>l/S ^ 



and 



T/= [\a^{t) [ d% [ d^w e-'{-*-+^*^)WA*.+v*A^) f dx ^,^(,,.,+,^) 

Jo Jc" Jc" J\x\<5 ^ 



vanish in the limit 6^0. 
For real r, p and q define 



^^(r,p,g) = / — e*^(P^+^) = / — e*-(«^+^). (5.10) 



L 



By integrating by parts, 

9L{t;p,q) = — — e*^(P^+^) + — / ^ dx + ^ ^ dx, 

iprL tpr Jj^ x"^ q X'^ 

and, therefore, for L > we have 
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Obviously, 

J\x\>l/6 X s ^ 



'\x\>l/5 

and, by dSHH), 



45 w*w5'^ 

< ^ + . 

v*vvt v*v 



'\x\<6 X 

Therefore 

45 w*w6^ 



As the function is locally integrable with respect to rf^ for n > 2, we conclude that 

J/ = 0(5) when 5^0. (5.12) 

Similarly 

J\x\<S X s « 

and repeating the above argument one obtains that 

J/ = 0(5) when 5^0, 

so that both X/ and X/' vanish in the limit 5 — > 0. This completes our proof of the first 
part of Theorem |21 

It is worth mentioning another formula for the regularized average of the inverse 
spectral determinant. 



R,{A,A*) = in-l) [ 



{1 - \z\^)''-^d^z 

l^|2<i det[e2j + (1 - z)I + (1 - z)AA* 



which is almost an immediate corollary of ()3.33|) and the representation 

1 1 /" ^-v''[e^I+{I^AU){I-AUr]v^2^ 



det[£2/„ + (/„ - Af/)(4 - AU)*] vr' 

1 



^-[v*{l+e'^)v+v*AA*v] ^tr{vv*AU+U''A*vv'')^2^ 



This formula, however, does not seem to be easy to handle in the limit e ^ 0. 

We now turn to the integral in ()5.8p and evaluate it in the limit e —>■ under the 
assumption that AA* has no repeated eigenvalues. The limit in ()2.15|1 follows immediately 
from the asymptotic relation ()5.2()|1 which is the end-product of our calculation. 

The following identities, which can be obtained from the Lagrange interpolation for- 
mula, see, e.g., [10], will be useful for our purposes. 



24 



Proposition 8 Suppose that Xi, . . . ,x„ are pairwise distinct. Then 

" 1 " 1 1 

n X - — ^ ~ ^ — f n T,. — T..' (5.13) 
j=i . 

and, for non-negative integer r, 



. , Xi — t 



where the hr, r = 0, 1, 2, . . . , are the complete symmetric functions 
It follows from ()5.13|) that 



det[AA* + - t)I - ieVi{x + i)/] 



where ai, . . . , a„ are the eigenvalues of AA* and 

1 

w{a,t,x) = -= 

a"^ + — t — ie\Jt{x + ^ 

By the calculus of residues, 

1 /•+~ dx 



w(aA,x) = — , (5.16) 

v/(a2 -t -£2)2 + 4^2^2' 



and putting ()5.15p and ()5.16|) into ()5.8|) we arrive at the following expression of Re{A, A*) 
in terms of the eigenvalues of AA*: 



n 

RM. A*) = J2 Fe{ai) n (5-17) 

where 



3 = 1 ky^J ^ 1 



Fe{a) = f ==, da^{t) = {n-l){l- tT-'dt. (5.18) 

This formula is convenient for finding Re{A, A*) in the limit e —>■ 0. 

If AA* > I, by letting e — > in ()5.17p and recalling ()5.13p we immediately obtain 

\imRM.A*) = f \ /r?^^^ 
£-0 ' ^ det(v4A* -t/)' 

thus reproducing the corresponding formula of Theorem^ part (ii). If A A* < / or if A A* 
has eigenvalues on each side of = 1, evaluation of the right-hand side in ()5.17|1 in the 
limit e ^ requires some work. 
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The integral in ()5.18|) is standard. There are different methods available to evaluate 
it. None seems to give an explicit expression for all parameter values. However, we are 
only interested in e —>■ 0, and in this regime 

Fe{a) = (n-l)(l- 02)^^-2 [7,_2Sgn(a2-l) + Lo(£,a)] +g„_2(a') + 0(£), (5.19) 

where 



In ^-f 


if a2 


<1, 


In 


if 


>1, 


In^ 

£ 


if 


= 1; 



Lo{e,a) 

7„_2 is the partial sum of the harmonic series, 

n . 



n-2 

7n-2 



sgn is the sign function, sgn(a;) takes value 1 if a; > 0, -1 if x < and if x = 0, and 
Q'n-2('2^) is a polynomial of degree n — 2 in with coefficients which do not depend on e. 
Details of derivation of ()5.19|) are given in Appendix ^ 

Let us now put into ()5.17|1 . In view of ()5.14j) the polynomial g„_2 gives no 

contribution. After rearranging the remaining terms we obtain, 

Rs{A, A*) = a{AA*) ln\ + p{AA*) + 0{e), (5.20) 



£2 



with the coefficients a and (3 given by 



1 

a{AA*) = (n-l)5^(l-a^2)"-^^(l-a^2) n ^ 

j=i k^j ^ 

n 

f5{AA*) = (n - 1) 5^ (1 - a|)"~^ V^(a?) n 
where 6 is Heaviside's step function 



j=l k^j '^fc '^■J' 




e{x) 



and 

{7„_2 + Ina^ - ln(a2 - 1) if > 1, 

-7„-2 + ln(l-a2) ifa2<l 

ln2 ifa2 = l. 



26 



As one would expect, the coefficient ot vanishes if A A* < / or AA* > I. This follows 
from identity ()5.14|) . If AA* > I then the constant 7n-2 gives no contribution, again by 
flQIjl and 



" 1 
limi?,M,A*) = (n-l)y (l-a2)"-2 In^^ n 9 , , ... 



don{t) 



j=i "i ~ ^ k^j "fc ~ "i -^0 

Similarly, if AA* < I then 

n 

Thus, fl5.8|l indeed reproduces formulas of Theorem [T] part (ii). 



dcTnit) 



6 Rank-one perturbations of CUE and GUE 

In this section we express the mean eigenvalue density for the random matrix ensembles 
()1.4|) and ()1.6p in terms of the spectral determinants. Our calculation is inspired by similar 
calculations in fr7\ ITHj and makes use of a process known as eigenvalue deflation which 
was introduced in the context of random matrices in |44j. 

We need to recall a few facts about elementary unitary Hermitian matrices Let 
V be a column- vector in C". The matrix 

Rv = In — 2vv* /\v\^, = v*v, 

where /„ is identity matrix, defines a linear transformation which is reflection across the 
hyperplane through the origin with normal It is straightforward to verify that 

is unitary and Hermitian, 

Rv = Rv and RvRv = R^ = In- 

In the context of numerical linear algebra the matrices R^ are known as Householder 
reflections. Any matrix can be brought to triangular form by a succession of Householder 
reflections. We only need the ffist step of this process which we now describe. 

Let Wn be an n X n matrix and z and ) be an eigenvalue and eigen- 



vector of Wni so that 



W„x = zx. 



Without loss of generality we may assume that Xi > and |a;p = x*x = 1. Let Ci = 
(1,0,..., 0)^ and 

x + ei x + ei 

V = = . . (6.1) 

|a; + ei| v^2(l + xi) 
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Since the vector v bisects the angle of x and ei, we have R^x = —e^ and -Rt,ei = —x. 
Therefore R^WnRvei = ze^ and (recall that R^ = In) 

Wn = R. (l ^* ] R., (6.2) 



n— 1 



for some Wn-i and w. Note that Wn-i is {n — 1) x [n — 1) and w* is 1 x (n — 1). 
Obviously, applying this procedure again (to the matrix Wn-i) and again, one can reduce 
Wn to triangular form by means of unitary transformations. Such factorization is known 
as Schur decomposition. 

It is convenient to write v = {vi,q)'^, where q = {v2, . . . ,Vn)'^ . Since i; is a unit 
vector, vl+ |gp = 1. Note that the first equation in ()6.1|) reads vi = a/(1 + xi)/2. Since 
< Xi < 1, we must have 1/2 < fi < 1. Therefore, 

= y/l - |q|2 and ^ < \q\^ < 1. (6.3) 
In terms of q the matrix R^ is given by 



= i'^ ; ,'1 , = ( i^^r'i . (6.4) 



2v^q In-i-2qq*J \-2y/l -\q\^q I^.i - 2qq* 

The incomplete Schur decomposition ()(i.2j) gives rise to a new coordinate system in the 
space of complex matrices, the new (complex) coordinates being z, w, q and the matrix 
entries of Wn-i- There are no restrictions on the range of variation of z, w and Wn-i, 
and, in view of ()6.3p . the vector q is restricted to the spherical segment | < |qrp < 1. The 
Jacobian of the transformation from (Wnjk) to this new system of coordinates^. 



j,k=l 



dzdz ^ dwjdwj j-r dqjdqj ^ d{Wn-i) jkd{Wn- 



i)jk 



^, . azaz T— I- uwjUWj -i— r aujauj -w-r 

J{z,w,q,Wn-i) ^ [[ 2 11 

is given by (cf. Lemma 3.2 in |TH] ) 

J(z, w, q, Wn-i) = 22«-2| det(^J„_i - Ty„-i)P(l - |qn«-2(2|g|2 - 1). (6.5) 

We derive ()6.5p in Appendix iBl 

Suppose that we have a probability distribution 



«J=1 



•^For Jacobian computations it is convenient to consider z and z as functionally independent variables, 
so that (Pz = dRe zdlmz = dzdz/2. 
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on the space of complex n x n matrices. Then, following the argument of ^8], see their 
Lemma 3.1, the mean eigenvalue density Pn{x, y), z = x + iy, of Wn is given by 

Pn{x,y)= j (fWn-i j d^w j d'^q J{z,w,q,Wn-i)p(^R^(j^ ^ ^ R^^ , (6.7) 
where 

Since we integrate in the q-space over the spherical segment 1/2 < |qrp < 1, it is 
convenient to introduce spherical coordinates, 

q = Vi(T, t = |qp, cr = q/\q\. 

The element of volume in the g-space is then 

d^q=^t''-^dtdS{a), 

where dS{(T) is the element of area of the sphere |crp = 1. The range of variation of t is 
1/2 < t < 1. Next, on making the substitution 



r = (2^-l)^ i2t - l)dt = \dr, (1 - t)t = 1(1 - r), ^ = 
the expression for the Jacobian becomes simpler, 

J{z,w,q,Wn-i)d'^q = J{z,w, Vi(T,Wn-i)d^q 

= 22"-3| det(;zJ„_i - Wn^i)\^[{l - t)t]"-'(2t - l)dtdS{a) 

= 1| det(^/„_i - Wn-i)\\l - ry-^drdS{a). 

Substituting this into ()6.7p . we arrive at the desired formula for the mean density of 
eigenvalues in the ensemble with matrix distribution ()6.6|) : 



1 



z w 



Pn{x,y) = - Jd^Wn-iJ d^wJdS{a)J{l-rr-^dr\ det(z/„_i-iy„_i)|^p \R (^^ w^ J^ 

C2{n-1)2 Cri-1 \cr\^ = l 

Here 

a/t ^/l — r (T* 

R={ , 1- (6.9) 

-V^r^cr - (1 + ^)crcr* 
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We shall now apply this result to express the mean density of eigenvalues in terms of 
the absolute square modulus of characteristic polynomials for two ensembles of random 
matrices. 

Rank-one deviations from unitarity. Let f/„ be an n x n unitary matrix and 



= ( /! J • " s ^ s 1- 

The Haar measure on U{n) induces a measure on the manifold W*Wn = in the space 
of complex nxn matrices via the correspondence Wn = UnGn- The corresponding matrix 
distribution is uniform and can be conveniently described via matrix delta-function 



dP{W^^) = cJ{W:W^ - Gl) n '^^^"^f^^"^'^ (6.11) 

For Hermitian H, we define 5{H) as 

5{H) = \[5{H,,)\[5{H,,)5{H-u). (6.12) 

The normalization constant c„ can be easily computed by changing to the matrix 'polar' 
coordinates, 

2" l!2!---(n-l)! 

Note that the eigenvalues of GnUn and UnGn coincide, and, therefore, for the purpose 
of calculating the eigenvalue statistics the ensembles GnUn and UnGn are equivalent. In 
view of ()6.2|) it is more convenient to work with matrices UnGn- 
Changing the coordinate system to 2, r, cr, Wn-i, 

\det{zIn-i-Wn-i)\'^d^z (1 -r)"~2rfr dS{(7) d^w d^Wn-i, 

where we have used the unitary invariance of the matrix delta-function. 
The matrix inside the delta-function in ()6.14j) is 



|zp — 1 + 7r zw* — 7a/ (1 — r)r cr* 

zw - 7a/(1 - r)r cr W*_iWn^i - I„-i + ww* + 7(1 - r)cr(T* 

and the delta-function factorizes into the product of the three delta-functions, correspond- 
ingly. On substituting this into ()6.8|) we obtain 

Pn{x,y) = ^ j\ det{zln-l - Wn-l)\'f{Wn-l) d'^Wn^l (6.15) 

C2(n-1)2 
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where 



J. 

fiWn-i)= jdS{a) j{l-rr-^drS{\z\^~l+jr)x 

|<t|2=1 

J £w 5 {W:_,Wn-i-In-i + ww*+^{l-r)(TCT*) 6 (^zw-^^/{l 



, ... ... .V , , y , V V r]r <T ] . 

C"-i 

Note that 

7a/ (1 — r)r cr 



Therefore the integral over w yields 
1 



|^|2(n— ^ ( W:-l^n-l - In-l + 'ij-, ^ 



7(1 - r)(7r + 


\z\ 


?) 






2 



and 



1 

/(W^n-i) = jdS{a) j{l-rr-^dr5{\z\^-l + ^r) 



X 



5 i^:_iW^„_i - /„-i + ^ — — ^ (TO- 



7(1 - r)(7r + 








\z\ 


2 



It is apparent that ifj-zp < 1 — 7 or \z\'^ > 1 then the integral over r vanishes and 
f{Wn-i) = 0. Therefore pn{x,y) = for such values of 2;. If 1 — 7 < < 1 then the 
integral over r produces a non-trivial contribution and 

1 r / ^ \ \'^\ / 1 I I 12 

|<t|2=1 

Introducing 

7 

we can rewrite the above expression in a shorter form, 

= U j J dS{(T) 6 {W:_,W^^r - 4_i + 70-0-*) . 

On substituting this into ()6.15|) . we arrive at 



7 - 


1 + 


\z\ 


2 




\z\ 


2 



n-2 

Cn / 7 
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1^2^ y c?2lV„_,|det(^J„_i-l^„_i)|2x 

y dS{fT) 6 {W:_,Wn-l - In-l + 1(7(7*) . 



C2(n-1)2 



|<tP=1 
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Since the matrix crcr* is unitary equivalent to the matrix 



1 

0„_i 

the integration over cr can be easily performed yielding 



Pn{x,y) 



2-f\z\ 



1 



det(z4_i - Wn-i)V5 Wl_^Wn-x - GiM'W, 



'n~l, 



C2(n-1)2 

where Gn-i is the (n — 1) x (n — 1) matrix (cf. (jfi.lOj) ^ 



Gn-l 



VT^ 



In-2 

and Vol(5'^"'^^) is the area of the unit sphere in M^("~^^ 



Vol(5 



2n-3\ 



It follows from (jfi.lHj) that Cn/cn-i = {n — l)!/7r". Hence 

CnVol(52"-3) n - 1 



(6.16) 



"Cn— 1 ; 



and finally 



Pn{x,y) 



n — 1 ( ^ 



n — 1 / 7 



n-2 



c„^i / |det(z/„_i-iy„_i)p5 hy^.^pv^.i-ctj SWn_, 



C2(n-1)2 



n-2 



det (2;/„_i - Un^iGn^i 



dUn~i 

^l\Zr \1J Juin-l) ' 

as claimed in ()1.5p . 

Rank-one deviations from Hermiticity. Let be a GUE„ matrix, i.e. random 
Hermitian matrix of size n x n with probability distribution 

and r„ be the n x n matrix 



j=l l<j<k<n 



^^^ = ^•0 oil' 
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Consider the random matrices 
Obviously, 

ReWr,:= ''^Y =Hn, and Imiy„ := " ^ / " = T,- 
The matrices Wn are complex and their probability distribution is given by 



d{Wn)ijd{Wn)ij 

Jl^iiii vvn — n) J_ J_ 

where C/3,„ is the normalization constant 



dP{W^) = c,,e-i-(^^^")^5(Imiy„ - T.) H ^ "^'f (6-17) 



%n=(- ^ (6.18) 



TT 



and (5 is the matrix delta-function ()6.12j) . 

Changing the coordinate system to z, w, r, cr and Wn-i, 

I det{zIn^i-Wn^i)\^dh (1 - r)"-2rfr c/^(o-) d'^w d'^Wn-i, 

where we have used the unitary invariance of the matrix delta-function. 
The matrix inside the delta-function in ()6.14p is 



Im 2; — 7r ^ + 2 a/(1 — r)r cr* 

^-^ + fv^(l-r)r a* Im - 7(1 - r)crcr*^ 

and the delta-function factorizes into the product of the three delta-functions correspond- 
ingly. On substituting this into ()6.8|) we obtain 

p4x,y) = \cp,ne-^ j I det(z/„_i - W^-i)? f{Wn-i)e--^'''^''''^"-'^' d^W^-i (6.20) 

C2{n-1)2 

where 

1 

\n-2 



/(iy„_i) = jdS{cT) J {l-r^-'dr 5{y-^r)5{lmWn-i-l{l-r)(T(T*)x 

|(t|2=1 



d^w e 4 6 { — + ^^{l-r)r a' 
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The integral over w yields j e ^'^'^^^ and we arrive at 



4 

1 



f{Wn-i) = \ JdS{(T)Jdr {l-rr~^e'''^'^'~'->6{y-^r)6{lmWn-i-l{l-r)(7(T'). 

|(t|2=1 

It is apparent that if y < or ?/ > 7 then the integral over r vanishes. Therefore 
Pn{x,y) = if 1/ < or y > 7. If < y < 7, then the integration over r produces the 

factor - and the constraint r = -, so that 

7 7' 



f{Wn-i) = J dS{a) 6 (Im W^-i - (7 - y)^^l ■ 

I,t|2_i 



On substituting the obtained expression for f(Wn-i) into ()6.20|1 . we arrive at 

_ ,,\n-2^-^-Pi'r-y)y 



Pn{x,y) = '''-^^ '~ I rfW„_i|det(z/„_i-Ty„_0|V§*^(^^('^-)^x 

C2{n-1)2 

J dS{cT) S{lmWn-l-{l~y)(TCT*). 

The integral over cr yields 

where r„_i is the [n — 1) x (n — 1) matrix 

and Vo^^*^""^) is the area of the unit sphere in R^("~^) ()6.16|) . We have that 



/ /3 \ / 



1/2 



and it now follows that 

^_ 1 - y)"-2e-¥-/^(7-y)y 

^' ~ 4^2^ (n- 2)! 



C/3,n-l 

C2(n-1)2 

as claimed in p.7j) . 



y d'Wn-i\ det(^J„_i - W^„_i)|V^^(^^(^"-)'(5 (im - f„^i) , 
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A Appendix 

In this appendix we evaluate the integral 

(1 - t)'' dt 



V(t-a2 + £2)2 + 4^2^2 



in the limit £ — > 0. 

We shall use the following fact from Calculus. If P{t) is a polynomial of degree k then 
(integrate by parts) 

J ^/t^+pt + q J + pt + q 

where Q is a polynomial of degree k — 1 and A is a constant. For Q and A one has the 
equation (differentiate ()A.1|) ) 

P{t) = Q'{t){e+pt + g) + lQ(t){e+pt + q)' + A. 
It follows from this that 

_n t=l 



Q.,a(t)V(t-a2+ £2)2+4^2^2 + A,,,Jo(£, a) 

J t=0 



= Q.,a(l)V(l-a2 + £2)2 + 4^2^2_ g^_^(0)(£2 + + A,,„To(£, a), 

and the equation for Q{t) and A is 

(1 - tf = g;,,(t) [{t -a' + ey + 4e'a'] + Q,,,(t)(t - + e') + A,,,. (A.2) 

It is apparent from (jA.2jl that Q{t) and A must be polynomials in and £2 and, therefore, 
in the limit e —>■ 0, 

QeAt) = Qa{t) + 0{e^) and A,,, = A, + ©(e^), 

and 

(1 - t)" = Q'^m - ay + Q^mt - + Aa. 

This equation for Qa and Af, can be explicitly solved, the solution being 

A. = (l-aY and Qa(t) = ^ (^) (t - a^y-^(l - a^"'. 

1=1 ^ ^ 

Note that at Qa{0) is a polynomial in of degree A; — 1, and 



1=1 ^ ^ 



-(l-aY-S. 
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where 7^ is the partial sum of the harmonic series, 

We now turn to 2o(£^, a^). RecaUing the table integral 

we have 
At = 1, 

For 0} 7^ 1, 



^ , 2 2^ , 1 - + + V(l - a2 + £2)2 + 4£2q2 

Jo(£ , a ) = In 



To(£', 1) = In = £ + 



^(l_a2 + £2)2 + 4^2^2 = |1 _ ^ ^'[^ ^ + 0(5^), 

1 — a 



and therefore 

l-a2 



2 „2\ _ ) £ 



In + 0(£2) if a2 < 1, 



In- + 0(£') ifa2>l. 

After collecting all relevant terms we arrive at the desired formula 

Xk{e\ a^) = (1 - a^fe{l - a^) In ^ + (3{a^) + qk{a^) + 0(£), 
where 6* is the Heaviside step-function, 

{1 if X > 0, 
I ifx = 0, 
if X < 0, 

qk is a polynomial of degree and 

(5{a^) = sgn(a2 - 1)(1 - a^fi^fk - In |1 - a^]) + ^(a^ - 1) Ina^. 
We use the convention according to which the sign function, sgn(a;), vanishes at 
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B Appendix 

In this appendix we derive equation ()fi.5|) . Let 



Wn-iJ ' V-Vl - <1*Q<1 In-i~2qq* 

When z,w,r] and Wn-i get infinitesimal increments dz, dw, dq and dWn-i the matrix 
Wn gets increment 

dWn = dR(f, ]r + r f "If j!!"* ]r + r(^^ ] dR. 

Since i? is unitary Hermitian, the matrix RdR is skew- Hermit ian, so that 

dT = RdR = ( / ,m 
for some /, /i and T„_i. Also = —{dR)R, and it follows that 
RidW^)R = dT(' ,f]dT+(i' 

Zl-Wn-I)dh dhw* + dTn-lWn-l-Wn^ldSn-l dWn-l^ ^^ 

Let dM = R{dWn)R- It is apparent that 

n n 

W d{Wn)jkd{Wn)jk = n dM.kdMjk. (B.2) 

j,k=l j,k=l 



On the other hand, it follows from (jB.lj) that 

n n— 1 n— 1 n— 1 

Y[ dMjkdMjk=\det{zI-Wn-i)\'^ dzd^ Y[dwjdwj Y[dhjdhj JJ d{Wn^i)jkd{Wn-i)jk. 

j,k=l j=l j=l j,k=l 

(B.3) 

To complete our derivation we now compute the Jacobian of the transformation from 
h to q. Recall that dh is the (2,l)-entry of the matrix dT = RdR. A straightforward 
computation yields 

dh = (2a + b){dq*)qq + a{dq) + bqq*{dq), (B.4) 

where 

1 - 2q*q 



a = — 2y 1 — q*q, b 



VI - q*q 

Equation ()B.4jl can be written as 

dh = {al + bqq*){dq) + (2a + b)qq^{dq), 
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and, therefore, 

^dq\ _ f al + bqq* (2a + b)qq'^\ f dq 



dq J \^(2a + h)qq* al + 6qq J \dq 
It now follows that 

n—l n—l 

dhjdhj = det(a/ + L) JJ^ dqjdqj (B.5) 
i=i i=i 
where L is the 2{n — 1) x 2{n — 1) matrix 

bqq* (2a + b)qq'^ 
(2a + b)qq* bqq^ 

If we find the eigenvalues of L, we shall know det(a/ + L). 

To solve the eigenvalue problem for L, we observe that if {f,g)^ is an eigenvector of 
L then 

b{q*f)q + {2a + b){q^g)q = Xf 
(2a + b){q*f)q + b{q^g)q = Xg 

for some A. If A 7^ we must have f = Ciq and g = C2q for some Ci and C2, and 

b{q*q)ci + {2a + b){q*q)c2 = Aci 
{2a + b){q*q)ci + b{q*q)c2 = Aca 

This reduced eigenvalue problem yields the two non-zero eigenvalues of L, 

Al = —2aq*q and A2 = 2(a + b)q*q. 

It is now apparent that A = is an eigenvalue of L of multiplicity 2{n — 2). This fact 
can be verified independently of the eigenvalue count by observing that for any vector u 
which is orthogonal to q, 

Ko)=°-"(^°)=°- 

It follows now that 

det(a/ + L)= a2("-2)(^ + Xi){a + A2) = i-2Y'''\l - q^qT'^il - 2q*q). (B.6) 
Collecting (|R2|l - (|R8|) and (jR5|) - one arrives at ^!B>- 
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